Spin-Orbital Fluctuations and a large mass enhancement in LiV204 
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We present a scenario that the muhi-component fluctuations, especiaUy those of the spin-orbital 
coupled modes, lead to the mass enhancement observed in LiV2 04. This phenomena is possible 
because all these modes are fluctuating due to the geometrical frustration. To illustrate this mecha- 
nism, the t29-orbital Hubbard model on the pyrochlore lattice is studied based on the random-phase 
approximation. We derive the generalized susceptibility in the SU(6) spin-orbital space and calcu- 
late the free energy by using a coupling-constant integration. The estimated speciflc heat coefflcient 
is of the correct order of magnitude to explain the experiment. 

PACS numbers: 71.27.-|-a, 71.28.+d, 74.20.Mn 
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I. INTRODUCTION 

The metallic spinel LiV2 04 is a unique Sd heavy- 
fermion (HF) compound [Q. Its 7 coefficient amounts 
to 420 mJ/mol-K^ and, strikingly, the low-temperature 
properties below T* =2Q ^ 30 are quite similar to those 
of the lanthanide or actinide HF compounds j^, |2|. Since 
the d electrons, 1.5 per vanadium ion, are responsible for 
both transport and magnetism in LiV204, it is not trivial 
whether one can apply a scenario analogous to the con- 
ventional HF mechanism. Moreover, the large spatial ex- 
tent of d orbitals compared with well-localized / orbitals 
is also unfavorable for the electrons to be treated as lo- 
calized. In fact, some transport properties, such as the 
steep increase of resistivity above T* g and the pressure- 
induced metal-insulator transition are different from 
those of the typical /-electron HF compounds. 

Apart from the HF behavior, geometrical frustration is 
another specific feature of this system. The B site of the 
normal-spinel (AB2O4) is known to form the network of 
corner-sharing tetrahedra (pyrochlore lattice). In most 
spinels, the geometric frustration is partially released by 
the structural phase transition. As a result, the ordered 
ground state is realized in the tetragonal phase, such as 
the possible charge ordering in Fe304 or the Neel order- 
ing in ZnV2O4 0. In contrast, LiV204 remains cubic and 
no static magnetic order was observed down to 0.02 K. 
The Curie- Weiss fit of the susceptibility (100 <T < 300) 
results in an effective spin- 1/2 moment per V with an 
antiferromagnetic (AF) coupling (0cw = ~37i4r)Q. AF 
long range order (LRO) will be suppressed by the geomet- 
ric frustration and, indeed, the inelastic neutron scatter- 
ing experiment captures the development of short-range 
AF correlations below 0cw||^- Interestingly enough, the 
large enhancement of 7 is shared by other geometrically 
frustrated systems with strong AF spin fluctuations, like 
Y(Sc)Mn2|§| and /3-Mn(7|. These facts seem to suggest 
that the formation of the heavy- mass d electron is related 
to the geometric frustration. In view of the geometric 
frustration, it should be noted that, already at the 10% 
electron doping, the HF ground state is easily destroyed 



into the spin-glass phasedj. 

In LiV2 04, the crystal filed around each V atom is 
nearly cubic with a slight trigonal distortion. The band- 
structure calculations indicate that the Fermi level lies 
within the cubic-i23 manifold with a bandwidth of about 
2 eV|^, |9[ |ll|, which is actually the superposition of 
the two components; the relatively narrow Aig (~ leV) 
and the broad Eg orbital (~ 2 eV) of the D^d group. Mat- 
suno etal. considered this trigonal split less important, 
and stressed the importance of the geometric frustration 
and/or orbital degeneracyl^. On the other hand, the 
LDA-I-C7 calculation by Anisimov etal. suggested that 
the correlation effect leads to the orbital polarization, 
where the Aig orbital is occupied singly to form a local- 
ized spin-1/2 moment and the rest 0.5 electron fills the 
Eg band. Thereby, they proposed the Kondo scenario 
based on the effective Anderson impurity model 

Followed by these band-structure calculations, various 
theoretical models are proposed in connection with the 
microscopic origin of the 'id HF behavior |lj, ^ |lj, |l^, 
[T^ , [l7[ . Especially, the role of the geometric frustra- 
tion has been featured in recent studies. Motivated by 
the Curie- Weiss law around room temperature, the for- 
mation of the localized s = 1/2 is assumed by several 
authors on the grounds of the orbital polarization ||l5[ |l^] 
or the local valence fluctuations[p7t. In these treatments 
the spin frustration is their main concern. 

In this paper, we consider that the present system re- 
alizes the prototypical itinerant frustrated model and in- 
vestigate the t2g-band pyrochlore Hubbard model with 
particular emphasis on orbital fluctuations. Within the 
random phase approximation (RPA) scheme, all kinds 
of fluctuations are included simultaneously without prej- 
udice. In principle, geometric frustration prevents any 
kinds of LRO and thus the large 7 would be expected 
by the resultant enhanced fluctuations. In this scenario, 
not only the spin but also both of the orbital and spin- 
orbital fluctuations are responsible for the enhancement. 
To our knowledge, this is the flrst theoretical attempt 
to attribute the heavy-mass d electron directly to the 
orbital degrees of freedom. The instability itself of the 



2 



same model was discussed by Tsunetsugu|18 . We have 
explicitly calculated the T-linear coefficient of the imagi- 
nary part of the generalized susceptibility and found the 
large enhancement of 7 mainly due to the spin-orbital 
fluctuations. 



II. MODEL HAMILTONIAN 
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The pyrochlore lattice is a fcc-array of tetrahedra. The 
cubic unit cell shown in Fig. |l] contains 16 lattice points 
covered by 4 tetrahedra and the primitive unit cell in- 
cludes a single tetrahedron. Therefore the t2g electrons 
are specified by the following 4 indices; the unit cell (de- 
noted by j = 1, • ■ • , 4L'^), the sublattice (n — 1, • • • , 4), 
the orbital {m = d^y, dyz, dzx), and the spin {a / I), 
where L is the number of the cubic unit cell along one 
direction. By using the standard notation of the multi- 
band Hubbard model, the t2g-band Hubbard model on 
the pyrochlore lattice is given hy H = Ho +Hi; 
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Here 



is the transfer integral between the m or- 
bital at n sublattice and the m! orbital at n' sublat- 
tice when (m,n) ^ [m',n'). In addition, the parame- 
ter of t™™ (= ±D/3) type is also included to describe 
the trigonal split of i2g, ^aig — = D. Since the 
nearest-neighbor (NN) hoppings already depend on the 
bond direction and the orbital symmetry the NN tight- 
binding model would be sufficient for the present study. 
By the symmetry relation, all the hopping matrices are 
represented by the following three independent parame- 



ters, to = h = tll^\ and ^2 = qf^- The Ta 



symmetry of the unit cell, as well as the spin SU(2), are 
incorporated in the 12 x 12 matrix Tio, though here we 
do not display the matrix explicitly. 

The tight-binding fit of the LAPW band-structure 
calculation gives the parameters in Tio to be to — 
-0.281 eV, ti = -t2 = 0.076 eV, and D = 0.140 
eV0. By applying the orthogonal transformation Cfc^CT = 
J2v Ufj.^{k)aki,a-, the single-electron Hamiltonian is diag- 

onalized into Ho = J^kau <^kua\vaO'kvc!, where ^ = (m,n) 
and u are 12-dimensional indices. The calculated t2g- 
multiplet (e/cjy, v = 1, • • • , 12) are shown in Fig. 3 along 
some typical symmetric lines. The almost flat bands near 
the Fermi level results in high density of state (DOS) 
at the Fermi level. Actually, the 7 coefficient estimated 
from the DOS of the band-structure calculation is rel- 
atively large; 7band ~ 17 mJ/mol-K^. However, the 
experimentally observed value (7cxp) is about 25 times 



FIG. 1: The cubic unit cell of the pyrochlore lattice (left 
panel) and the three independent hoppings (right) 
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FIG. 2: Band structure of the t2g multiplet for the present 
parameter along symmetric lines of the fee Brillouin zone. 
The horizontal straight line is the Fermi level for the quarter 
filling. 



larger than 7band. As a possible origin of this enhance- 
ment, in the following, we consider the effect of electron 
correlations from the weak-coupling limit. 



III. FREE ENERGY AND RPA EQUATIONS 

In order to describe various kinds of fluctuations con- 
cerning the t2g electrons in a simple way, we introduce 
the thirty- five generators of the SU(6) group X^{'^ = 
1, • • • , 35) and the identity operator X'^, where the nor- 
mahzation condition is Tr X'^X'<' = 6''''' /2. The SU(6) 
generators are classified into the pure spin, the pure or- 
bital and the spin-orbital coupled components. Spin and 
orbital degrees of freedom are described by the gener- 
ators of SU(2) (£""; three dimensions) and SU(3) (T'^; 
eight dimensions) groups, respectively, and the spin- 
orbital coupled modes are made of the products of the 
two {2^/3S"T'^; twenty four dimensions). To be explicit, 
these operators are represented by using the electron cre- 
ation and annihilation operators as follows. 
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TABLE I: The matrix elements of the diagonal interaction 



matrices, A' 



(c/s) 
/3/3' 



for each irreducible representation. 



Ai(/3 = 0) ^(1,4,6) ^(2,5,7) 1-3(3,8) 



-U - AU' + 2J 
(7 + 2J 



?7' - 2J 

U' 



U' - 2J 



-U + 2U' - J 
U-J 



where a and A are the PauU and the Gell-Mann ma- 
trices in the standard notations jl^, respectively. Note 
that the overall signs of non-diagonal operators for dif- 
ferent sublattices are selected to be consistent with the 
Td unit-cell symmetry. Since the trigonal distortion in 
Tio violates the 0(3) orbital symmetry, the usual repre- 
sentations of the t2g orbital by the fictitious / = 1 and 
its higher moments do not reflect the proper symmetry. 
Instead we construct the tensor products of the real I — 2 
moment and reduce them into the irreducible represen- 
tations of Td in the cubic t2g subspace. As a result, we 
find that A is the proper set of basis, which is classi- 
fied into three irreducible representations; I^g-dipolar mo- 
ments (A2, A5, Ay), Pg-quadrupolar moments (Ai, A4, Ag), 
and /-g-quadrupolar moments (A3, Ag). 

Next we transform the interaction part of Hamiltonian 
in terms of X~''s so that 
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except for a constant energy shift. Here a-y^'s are a func- 
tion of [/, U' and J. Such transformation is always pos- 
sible for any U, U' and J, though the choice of a^^ is not 
unique because of the spin rotational invariance and/or 
the identity = njmna- By applying a couphng con- 

stant integration together with Eq. (|^), the Free energy 
for this system is given by 

F(U,U',J) = F{Q) + — / dujcoih'— 
27r y_oo 2 
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where the momentum sum is taken all over the 1st Bril- 
louin zone and the generalized susceptibility per tetrahe- 
dron is defined by 
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dre^(-+^^)-([X,^„(r),XZ;„, (0)]\.(9) 



We apply the random phase approximation (RPA) in 
the evaluation of the generalized susceptibility (x). In- 
cluding the charge degree of freedom (X°), our x is the 
144 X 144 matrix with the row index nj and the col- 
umn n'y . Since we are concerned with the paramagnetic 
state, as suggested by the experiment, the SU(2) symme- 
try ensures that x is decomposed into the spin singlet and 
triplet sectors. Each sector is represented by the charge- 



charge (denoted by x' 



1313' (c 



(x^f,^*''), respectively. In these two 36 x 36 matrices, the 
spin indices (a) is reduced and the rest indices are the 
row index nj3 and the column n'/3'. Here /? = is de- 
fined to represent the totally symmetric orbital operator. 
To sum up, Xnn'^'^^ describes the correlations among the 

charge (/3 = 0) and the orbital (/3 > 1), while Xnn'^^^ 
among the spin (/3 — 0) and the spin-orbital (/3 > 1). 
By using these notations, the RPA equations are given 
as follows. 



PP'{c/s) 



Ann' 



E 



/3/3i(0)r,A(c/s) /31/3'(C/S) . „^ 



Since the interaction is local, the interaction matrices 
(^(c/s)') g^pg always n-diagonal and independent of n. 
Thus the rest indices are row f3 and column /?'. Moreover 
the 9x9 matrices A^'^/") may have a simple representa- 
tion. This is because 7io has the Td symmetry and Eq. 
( pO| ) is invariant under any Td transformation. Therefore, 
AT^/s) 's must be diagonal and these matrix elements are 
the same within each irreducible representation, see Ta- 
ble |. 

Xnn'''^^ in Eq. ( |lo| ) is the paramagnetic susceptibility of 
the Hartree-Fock approximation. When we calculate x''"'' 
of the usual single band Hubbard model, the Hartree- 
Fock terms just renormalize the chemical potential. In 
our model the Hartree term similarly changes the diag- 
onal part of Tio, however the Fock terms introduce off- 
diagonal terms of the orbital hoppings on the same site, 
giving rise to the crystal field effect. This crystal field is 
totally symmetric concerning the entire unit cell, but its 
local symmetry on each sublattice is trigonal. The effect 
of this additional field is included in the band-structure 
calculations. Therefore the free paramagnetic suscepti- 
bility calculated by using the tight-binding parameters 

gives the appropriate Xnn'^^^ ■ actual calculation, 

we used the free x^^^ P^r tetrahedron of the form 



y^°^ (a to) 



4L3 ^ ^ 



k 1^11^2 



and the S^-S^ susceptibilities 



ikut - f-k+qi^2 + K'^ + 

X {k)U^^^^ {k)U^^^^ {k -f q)U^^^,^ {k + q), (11) 

and rotated the orbital basis from {mim2m^mi\ into 
{/?/?'}, where /(e) is the Fermi distribution function. 
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IV. ESTIMATION OF THE SPECIFIC HEAT 
COEFFICIENT 

Since we are primarily interested in the mass enhance- 
ment in 7, let us now focus on the w-linear part of Imx 
in the Free energy of Eq. (||) . As is seen from the band 
structure (Fig. y), nothing singular would happen in 
X^^), therefore a constant i?ex(°) and an a;-linear /tox'") 
behaviors are expected in the leading order of co. Then it 
is sufficient to solve the RPA equations in the lowest or- 
der of /tox'"' oc iu. For the 1/2-filled s-band pyrochlore 
Hubbard model, in contrast, the complete matching of 
the Fermi level and flat bands brings a non-Fermi liquid 
behavior pC|, |2l[ |. As for the still- unfixed a^^, we take 

- -^A^^ = -i (a(^) e a(^) e a(*) e a(^)) . (12) 

It is easy to check that the a^-y's of this form reproduce 
the original interaction term, Eq. (0). This choice is 
natural in the sense that A-y-y reflects the symmetry of 
the non-interaction term (SU(2) and T^). We neglect the 
weak dependence of x^^^ on the interaction through the 
Fock terms for simplicity, then the a;-integral in Eq. (||) 
can be performed analytically. After some calculations, 
7 coefficient {—~d^F / d^T) per tetrahedron is given by 

7RPA = E (^''^^^ + 37''(^)) , (13) 

with 

qn mpi 

x{l~-2K^;l;lRexi\T\<i.^)y' ■ (14) 

This equation is exact after taking u!^S,T and L 
oo limit. Note that we neglect the zero-point fluctuation 
term in the uj integral because it will only show a weak 
temperature dependence through that of x^"^ . 

Finally we numerically estimate 7rpa in the SU(6) 
limit (U = U' and J = 0), where the orbital fluctu- 
ations are mostly enhanced. For this purpose, we cal- 
culated the free susceptibility, Eq. (|ll|), for the pa- 
rameter set of {uj,S,T) = (0.02eV,0.02eV,0.02eV) and 
L = 16. We checked the convergence of x''"'' against lj 
and L along some symmetric lines. Especially, we con- 
firmed that Imx^^^ shows a clear w-linear dependence 
in the range of w = 0.01,0.02,0.04,0.08 eV. The con- 
vergence of the cj-linear coefficient regarding the system 
size {L = 12, 16, 24, 32 ) is fairly good except for a few 
points very close to the F point. As a result, we ob- 
tained the 7-value per LiV204 mole as shown in Fig. 
|. The total 7 is the summation of the contributions 
from charge (denoted by 7c = 7°''^'), spin (7s — 37"^"'), 
orbital (70 — Y!ffi=i^^^^^)^ ^^'^ spin-orbital fluctuations 
(70 3X;Li7''^''^)- The critical value Uc is about 0.92 



eV and there the complicated eigenvectors, made of the 
four-sublattice summations of the spin and spin-orbital 
components with quadrupole orbital modes, characterize 
the instability at the wave vector near F point. Since 
wc fixed uj finite and took the limit of g — *■ 0, our calcu- 
lated i?ex^*''' does not correspond to the uj limit at 
F point. Therefore we can not determine that the most 
prominent modes are just at the F point or not. The real 
part of x^"' is studied for different hopping parameters 
at zero temperature jl^. In the present calculation, we 
do not detect the peak structure around X and L2 points 
found in ||l^ . The reason for the difference might be that 
our calculations have done at finite F and/or these sin- 
gle modes are quite sensitive to the nesting of the Fermi 
surface. 

Fig. m shows that the correct order of 7 can be repro- 
duced even away from the critical point, say ~ O.SUc- 
Around there, the enhancement of 7rpa (about 80%) 
is mainly attributed to the various spin-orbital fluctua- 
tions. Both spin and spin-orbital fluctuations contribute 
7RPA by the same order as a single component. How- 
ever, the magnitude of the components, 3 for spin and 
24 for spin-orbital modes, results in the different contri- 
bution in 7RPA- In order to gain the same enhancement 
solely by the spin fluctuations, U must be very close to 
Uc even after the inclusion of J term. Away from the Uc, 
the enhancement of 7rpa comes from the entire momen- 
tum space because the spin and spin-orbital fluctuations 
have rather weak g-dependences. This result seems to be 
consistent with the naive picture that the geometric frus- 
tration induces local fluctuations, which is responsible for 
the heavy mass. 
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FIG. 3: The graph shows 7 value per formula unit of LiV204 
mole. The contributions from charge, spin, orbital, and spin- 
orbital fluctuations are labeled by 7c,7s,7o and 7so, respec- 
tively. 
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V. CONCLUSION 

In conclusion we have studied the origin of the large 
enhancement of 7 in LiV2 04 based on the 1/4 filled t2g- 
band pyrochlore Hubbard model. We introduced T^i ir- 
reducible operators for orbital density operators and de- 
rived the RPA equations in a simple form. The 7rpa 
obtained by the coupling constant integration of Imx is 
enhanced typically by order of 10 due to the spin-orbital 
fluctuations compared with the 7s due to the spin alone. 
This conclusion itself seems to be general for the orbital 
disordered system and thus it may be possible to apply 
the present scenario to other systems, such as Y(Sc)Mn2 
and /3-Mn. To be more specific to LiV204, the calcu- 
lated 7rpa value is of the same order as 7cxp- How- 
ever, it should be noted that our simple theory does not 
include the effect of the geometrical frustrations suffi- 
ciently, since the mode-mode coupling between various 
fluctuations is neglected. In general, the RPA approxi- 
mation overestimates the fluctuations, therefore the most 
prominent mode itself would be suppressed and the Uc 
becomes larger if we include the couplings between the 
modes with different momenta. At the same time, the 



overall structure of % in g space would change into more 
and more structure-less and thus the additional enhance- 
ment of 7 could be expected. 

It is well-known that the mode-mode coupling theory 
leads to the Curie- Weiss behavior of the magnetic sus- 
ceptibility. As for the it is of order of 0.1 in the 
RPA approximation. This value would be increased by 
including J, although it is not clear whether the Hund's 
rule coupling is sufficient to obtain i?-vv ~ 1-7. To address 
the similarity to the 4/ HF compounds, we need to de- 
velop a more sophisticated theory including the coupling 
between fluctuations. The fluctuation exchange approxi- 
mation for the multi-band model would be useful as the 
next step. 
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